Non-Fermi liquid behavior in transport through Co doped Au chains 
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We calculate the conductance as a function of temperature G{T) through Au monoatomic chains 
containing one Co atom as a magnetic impurity, and connected to two conducting leads with a 
4-fold symmetry axis. Using the information derived from ab initio calculations, we construct an 
effective model Heft that hybridizes a 3d^ quadruplet at the Co site with two 3d** triplets through 
the hopping of ^Axz and 5dyz electrons of Au. The quadruplet is split by spin anisotropy due to 
spin-orbit coupling. Solving Hes with the numerical renormalization group (NRG) we find that at 
low temperatures G{T) — a — h\/T and the ground state impurity entropy is ln(2)/2, a behavior 
similar to the two-channel Kondo model. Stretching the chain leads to a non Kondo phase, with 
the physics of the underscreened Kondo model at the quantum critical point. 

PACS numbers: 73.23.-b, 71.10.Hf, 75.20.Hr 



Transport properties in nanoscopic systems are be- 
ing extensively studied due to potential applications and 
their role in basic research. In particular, the Kondo 
effect, in which the spin 1/2 of a molecule or a quan- 
tum dot (QD) is screened by one channel of conduc- 
tion electrons below a characteristic temperature Tr-, 
has been experimentally observed and the conductance 
as a function of ternperature Gi^T^ is in excellent agree- 
ment with theory [l|, l^]- In spite of the large effect 
of correlations, these systems are Fermi liquids where, 
for example, the conductance has the expected behavior 
G{T) ~ G(0) - aT"^ for T <^Tk with Tk the Kondo tem- 
perature and a some constant. More recently, the under- 
screenedKondo (UK) effect, leading to singular Fevini liq- 
uid physics [3|, y] has been observed, and quantum phase 
transitions (QPTs) involving partially Kondo screened 
spin-1 molecular states were induced by externally con- 
trolled parameters [5|-l7|. 

The overscreened Kondo effect, the simplest manifes- 
tation of which is the two-channel Kondo (2CK) model, 
is even more interesting because inelastic scattering per- 
sists even at vanishing temperatures and excitation ener- - 
gics. Therefore the system is a non-Fermi liquid, exhibit- ■ 
ing fascinating low-energy properties |8l4lC|. In particu- ■ 
lar, the impurity contribution to the entropy is ln(2)/2 • 
and the conductance per channel at low T has the form ' 
G{T) ~ Go/2±aVT, where Go is the conductance at zero ' 
temperature in the one-channel case. However, experi- ' 
mental observations of 2CK physics have been elusive, ' 
particularly due to its instability against the asymmetry ' 
of the channels [10|. For example, while 2CK physics is - 
expected to take place at the QPT in the model of two ■ 
spin-1/2 Kondo impurities [9|, |lO|, inter-channel charge. 
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Figure 1: (Color online) Structure of the systems studied, 
a) Au chain doped with Co. b) trimer with one Co atom 
connected to BCC leads. The shaded region indicates the 
unit cell used. 



transfer spoils the quantum critical point (QCP) [iC 
its observation in double QD systems |lll |. 
Oreg and Goldhaber-Gordon 



and 



12l | proposed to mod- 



ify the simplest setup of the one-electron transistor (a 
QD between two conducting leads) by adding a second 
large QD with level spacing (5 < Tk , which acts as a sec- 
ond channel of conduction electrons, with gate voltage 
fine tuned towards the QCP. Following this proposal, the 
differential conductance as a function of bias voltage Vf, 
was measured and a \/Vb behavior, characteristic of 2CK 



physics was observed for eVf, > fc^T [13 1. 

In this Letter, we show that 2CK physics is expected in 
transport through Co doped Au chains connected to con- 
ducting leads with a four-fold axis, as shown in Fig. \Tjp. 
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Figure 2: (Color online) Spin-polarized spectral density of '' 
different orbitals in a) the linear chain AuioCo depicted inioo 
Fig. [iK, and b) the system of Fig. [TJj. loi 



calculations for a Pt system with this geometry indicate 
that the x^ — y^ and xy orbitals are split by an energy 
Axy ss 1.4 eV, while the remaining orbitals arc unaffected 
[17| (for bulk gold, one has Axy ~ 1.0 cV). As shown 
in Fig. [2J3, we obtain a similar splitting here. We have 
taken the Au-Au and Au-Co distances as 4.93 and 4.53 
a.u, respectively. Choosing ep = 0, we obtain for the 
average energies of the different Co 3d orbitals E^^^y^ = 
-1.0 eV, E^^2_^2 = -0.1 eV (not shown), E^y = 0.2 eV 
and Exz = Eyz = 0.3 eV. The occupations indicate some 
admixture of the 3d^ quadruplets with 3d* triplets as 
described below. 

An important ingredient of H^s is the splitting D be- 
tween Co spin states with projection Sz = ±3/2 from 
those with Sz — ±1/2 originated from spin-orbit cou- 
pling (SOC), calculated as described below. 

The calculations above lead us to the following effec- 
tive Anderson Hamiltonian to describe the low-energy 
physics. 



In this case, the mobile electrons of the two channels are 
the 5d xz and yz electrons of Au. An advantage over 
the proposal of Ref. 



12 is that both channels are re- 



lated by SU(2) symmetry in the effective model Hcs and 
therefore fine tuning is not necessary. Furthermore, by 
stretching the chain, it appears possible to pass through a 
QCP with underscreened Kondo physics to another phase 
with reduced conductance, without Kondo behavior. To 
the best of our knowledge, our model has no analog to 
previous proposals for observation of 2CK physics [ij]. 

The strategy of this Letter is twofold. First we per- 
form ab initio calculations [15| to obtain an effective low- 
energy Hamiltonian. This is followed by a detailed analy-"' 
sis of the latter using analytical arguments together withios 
a numerical analysis based on the quasi exact Numericahog 
Rcnormalization Group (NRG) fl8M20|. no 

Spin-polarized spectral densities of different orbitalsm 
are shown in Fig. [2l For the linear chain (Fig. [2^), weiu 
find that Co is in a 3d^ configuration with total spin 3/2. lu 
Two 3d holes reside in the xz and yz orbitals. The re-ii4 
maining one is shared between the x^—y"^ and xy orbitals. ns 
The band dispersion (not shown) indicates that these or-ne 
bitals arc strongly localized, resulting in no Au orbitalsn? 



with 



y and xy symmetry at the Fermi energy ep- 



We obtain that the top of the 5d xz and yz bands ofng 
Au are slightly above ep (see the inset of Fig. [2^). The^o 
same result was obtained by previous ab initio calcula-121 
tions of pure Au chains 



23l |. While transport ex- 122 



76 pcrimcnts in Au chains indicate that these bands do not: 

77 cross the Fermi level [24|, it has been found that doping; 

78 with oxygen pushes these bands up establishing conduc-125 

79 tion through the ?)dxz 



and bdyz electrons of Au 



25,2 



80 Application of a gate voltage can have a similar effect. 127 

81 The degeneracy between x"^ — y^ and xy Co orbitals isi28 

82 broken in a system as represented in Fig. [T}). Previousi29 



i^eff = ^(£^3 + f M|)|M3)(M3| + ^2 51 l«A^2)(aM2| 
Mi a Ah 

+ E J2^.mM2<7\lM,){cl,^JaM2){M,\+R.c.) 



Here £"„, Mn are the energies and spin projections along 
the chain (2-direction) of the states with n E {2, 3} holes 
in the 3d shell of Co. The charge-transfer energy is de- 
noted by £^32 ~ E2 — Ej,. For n = 3, the state of maxi- 



7t 



xz\, yzX- X' 



^JO), where |0) 



mum spin projection is |4) 

is the full 3d^° shell and dl^ creates a hole with symme- 
try /? and spin a. Similarly, for two holes |al) = rfailf) 
{a <S {xz,yz}). The remaining relevant states at the Co 
site can be constructed using the spin lowering operator. 
The hopping amplitudes V^ between the Co atom and the 
states cli^^cr ^^^ assumed independent of fc, where cj^j.^,^ 
creates a hole in the 5d band of Au states with symmetry 
a at the left {i/ = L) or right (1/ = R) of the Co site. The 
hybridization of the Co atom to the mobile electrons at 
either side [23| is given by F^ = 2T:J2k l^i^l'^^i^ ~ ^uk), 
neglecting the small dependence on w. 

The Clebsch-Gordan coefficients (J2^M2cr| J3M3) lead 



32 



to a spin SU(2) symmetric Hamiltonian for I? = 
For any D, the model exhibits channel SU(2) symmetry 
28| which was also exploited in the NRG calculations 
below in. 

Depending on the value of D, the model exhibits rich 
physics. For D ~ 0, the model is an SU(2) underscreened 
impurity Anderson model (two channels with spin 1/2 
cannot completely screen a spin 3/2). For one channel, 
this model has been solved exactly [32|. In the Kondo 
limit (E32 ^ F) it displays singular Fermi liquid behav- 
ior [si y] observed in the conductance through molecules 
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containing Co"^ ions [6|. For two channels, similar physics 
is expected. For sufficiently negative D, \D\ ^ Tk, spin- 
ffip processes are inhibited, the Co spin projection S^ 
fluctuates between 3/2 and 1 or between -3/2 and -1. 
The properties of the system are described by two inde- 
pendent resonant levels [33|. For large positive D ^ Tk, 
we find that a Schrieff"er- Wolff" transformation leads to 
an anisotropic 2CK model for the interaction between an 
effective spin 1/2 for the Co states with Sz = ±1/2 and 
the spin s^ of the mobile states that hybridize with the 
Co at the Fermi energy 



H' 



2CK 



= jJ2(^^ 



■ 2SySy), 



(2) 



where J ~ ^{Vl + V^)/[E32 - (D/S)]. The critical be- 
havior of this model is the same as for the isotropic 2CK 
model 13J]. As a consequence, the system displays non- 
Fermi liquid behavior due to overscreening of the impu- 
rity spin. 

To calculate D wc exactly solved the 120 x 120 matrix 
of the Hamiltonian of the 3d^ configuration, as described" 
in Ref. [35|, including the effect of SOC via -ffsoc =" 
A^j Ij • Sj. To gain insight on the effect of the splitting 
Axy between 3d x^ — y^ and xy orbitals, we take the" 
above mentioned ab initio results for the energies of the 
system of Fig. [T]3, but assume E^y ~ —Ex2_y2 = Axy/2 
li^. For Axy = 0, we obtain D = -0.0426 eV, in-" 
dicating that the Co spin is oriented along the chain. 
As discussed above, no particularly interesting physics 
is expected in this case. However, as Axy increases, D 
also increases and changes sign for A 



xy 



0.8 eV. For 



Axy = 1.2 eV (corresponding to the setup of Fig. [TJd) we"^ 
obtain D — 0.0017 eV. One expects that Axy is reduced"" 
by stretching the device and a QPT would take place if 
£) = is reached 3l| . 

The total resonant level width F = F^ + r^ = 0.60 eV,"' 
is determined from the width of the peak of the xz 
and yz states above the Fermi energy in the ab initio 
calculations. From the position of this peak, we take 
£-32 = 0.3 eV, choosing the origin of one-particle energies 
at the Fermi energy ep ~ 0. For the numerical calcu- 
lations we take a band extending from ~W to W with 
W = 5eV. Repeating the ab initio calculations for the 
system of Fig. [T|3 either without the Au atoms in the^°'^ 
chain or without the left lead, we obtain F/j = 0.25 cV, 
rV = 0.34 eV, leading to a factor A = ATlTr/T^ = 0.97f °' 
[36| for the conductance [see Eq. (j3|)]. 

The effective model in Eq. ^ was simulated us- 
ing NRG, with the results presented in Figs. [3] and 2] 
[discretization parameter A and truncation energy Ek 
(in units of rescaled energies) are specified in the cap- 
tions]. In Fig. [3] we show the impurity contribution tO209 
the entropy. Si, averaging even and odd iterations. At2io 
high temperatures, the quadruplet and both triplets are2ii 
equally populated and 5^ ~ In(lO). As the temperature2i2 



w 




Figure 3: (Color online) Impurity contribution to the entropy 
as a function of temperature for D = (full line black), D = 
-0.0017 eV (full line gray) and D = 0.0017 eV (dashed line), 
calculated using NRG [A — 3, Ek ~ 6.8, keeping up to 6 200 
multiplets (15 800 states)]. 



is lowered below the charge-transfer energy E32 , the sys- 
tem is in the local moment regime characterized by a local 
spin 3/2, and therefore a plateau appears with Si ~ ln(4). 
For D = (black solid line), as the temperature is low- 
ered below the Kondo temperature Tk/W ~ 5.2 • 10~^, 
a partial screening of the local spin takes place and the 
ground state is a doublet with Si = ln(2). The same 
zero-temperature value is reached for negative D. 

In contrast, for positive D > Tk, Si starts to fall more 
rapidly for T < D, and at low temperatures it satu- 
rates at the peculiar value ln(2)/2 characteristic for 2CK 
physics. We have confirmed (not shown) that the same 
low-temperature value is obtained for D < Tk. From an 
analysis of the low-energy spectrum, we obtain that for 
Z? <C Tk the crossover occurs at an energy scale T* <C T*, 
which roughly goes as T* « TKexp[-(6TK/£')^/^]. A 
similar crossover scale has been observed in a one-channel 
spin-1 UK model with anisotropy [33] • The crossover 
scale T* is present for any finite D > 0. 

Non-Fermi liquid behavior is also encountered in the 
spectral density of localized electrons (not shown), which 
at zero temperature and at the Fermi energy in the 
Kondo limit is half the value expected from the Friedel- 
Langreth sum rule [10| . These peculiarities of the spec- 
tral density manifest themselves in measurable proper- 
ties, such as the total conductance through the system 
27] 



G(T) 



Ae 
~h 



dunrp{u:){-^), (3) 



with f{uj) — l/(e"/'"'^-^ + 1) the Fermi function, and 
p{u>) = J2aa P'^o'i^) *^^ impurity spectral density of 
states. 

The results for G(T) are presented in Fig. 01 For D — 







Figure 4: (Color online) Conductance as a function of temper-2( 
ature for D = (full black line), D = -0.0017 eV (full gray 2, 
line) and D = 0.0017 eV (dashed line), obtained with NRG^^ 
[A = 2, Ek = 6.5, keeping up to 12 300 multiplets (31 800 ' 
states)!. The inset shows a fit of the D > NRG data with 



G = a- h./T/T* with T'/W ^Q- 10"'' . The resulting fit is 
also replicated in the main panel (lower thick light line). ^'^ 



the conductance at zero temperature is G'(O) ~ 1.8 Go, ^™ 
where Gq = 2Ae^ /h. This is nearly the ideal value^''' 
2Go expected for perfect transmission with two channels."" 
However, with increasing temperature the conductance"^ 
falls more rapidly than for a usual Kondo model, reflect-"^ 
ing the physics of the UK model [3|47| . We have obtained"^ 
an excellent fit of the NRG data below Tk, using a phe-"" 
nomenological expression proposed for one-channel UK"' 
models ^ i3Sj (see fit Parks 2010 in Fig. El). 

As the anisotropy D is turned on, the conductance 
G{T) deviates from the behavior for D = for temper- 
atures T < \D\. In particular, for positive D the con-^^^ 
ductance at zero temperature, G(0), is nearly half the277 
corresponding value for D = 0. However, the most strik-278 
ing feature of G{T) for Z) > is that it decreases with"' 
increasing temperature following a \/T behavior, instead^"" 
of the characteristic T^ dependence of a Fermi liquid (as 
for D <). Fitting the results for T < 0.8 T* = 0.17 K, we^[^ 
obtain G/Go ~ 0.96 - 0A7^/T[K] (see inset of Fig. El). '^4 

The results above demonstrate measurable non-Fermi-^ss 
liquid properties in the model. We have checked that^"'' 
the results are robust and do not change if for exam-^**^ 
pie anisotropy in the triplets or more configurations of 
3d Co states are added. In real systems, an indepen-^g^ 
dent conduction channel involving 6s states of Au andjgi 
4s and Sd^z^-^^ states of Co contributes to the conduc-292 
tance, but this contribution is comparatively small be-^" 
cause it does not lead to a resonance near the Fermi en-^'" 
ergy (in contrast to the Kondo effect) and more impor-^*"^ 

296 

tantly, it does not alter the non-Fermi-liquid behavior. 
The latter is only affected by a magnetic field or physi-^gg 
cal ingredients that break the channel SU(2) symmetry,299 
as a splitting A^ between Co Sd^z and 3dyz orbitals.300 



This leads to a Fermi-liquid behavior below a tempera- 
ture T* « {vAz^Tk, where v = 1/2W w 0.1/eV is the 
density of states at the Fermi energy [39| . For our param- 
eters, Tk a; 2.6-10-^eV. Therefore, for moderate A^, T^ 
is very small and non-Fermi-liquid physics is observable 
in a wide range of temperatures. 

In summary, we have shown that it is possible to ob- 
serve non-Fermi-liquid behavior in transport through a 
Co doped Au monoatomic chain, in a setup that is not 
far from currently accessible technologies. It has the ad- 
vantage over previous proposals in similar nanoscopic 
systems that both channels are related by symmetry 
and therefore, they are truly equivalent. Moderate 
symmetry-breaking perturbations still allow to observe 
non- Fermi liquid physics in a wide range of temperatures. 
Furthermore, by stretching the system it appears feasi- 
ble to induce a quantum phase transition to a phase with 
singular Fermi liquid behavior and increased conductance 
(D = 0) and to another phase with reduced conductance 
{D<0). 
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